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We consider the condensate wavefunction of a rapidly rotating two-component Bose gas with an 
equal number of particles in each component. If the interactions between like and unlike species are 
very similar (as occurs for two hyperfine states of *^Rb or ^'^Na) we find that the two components 
contain identical rectangular vortex lattices, where the unit cell has an aspect ratio of \/3, and one 
lattice is displaced to the center of the unit cell of the other. Our results are based on an exact 
evaluation of the vortex lattice energy in the large angular momentum (or quantum Hall) regime. 



Experiments on rotating Bose gases have progressed 
rapidly in the last two years. Soon after the pioneer work 
at JILA and ENS [||, the MIT group created a vortex 
lattice with as many as 160 vortices |^. Recently, the 
JILA group has invented an ingenious method to increase 
the angular momentum of a condensate by performing 
evaporative cooling on a rotating normal cloud Q. In 
this process, the system spins faster and faster as it is 
cooled, while remaining close to equilibrium. With such 
rapid progress, one expects that equilibrium Bose gases 
with even larger angular momenta may be produced in 
the near future. 

At present, most experiments on vortex lattices are 
performed in single component Bose systems. It is nat- 
ural to ask what happens in two-component Bose gases, 
such as those made up of two hyperfine spin states of 
the same atom. The vortex lattices in such systems are 
bound to be more intricate than those in single compo- 
nent condensates, as the vortices in different components 
can move relative to one another to minimize the en- 
ergy. The purpose of this paper is to study the vortex 
lattices of two-component systems with large number of 
vortices, in what we call the "mean-field quantum Hall 
regime" . This is the regime where mean field theory re- 
mains valid so that each component (labeled by an index 
i = 1,2) is characterized by a condensate wavefunc- 
tion 4"^; yet the angular momentum of the system is so 
high that V&i is made up of the orbitals in the lowest 
Landau level in the plane perpendicular to the rotation 
axis. It has been shown recently that this regime will 
emerge in a three dimensional Bose gas at sufficiently 
high angular momenta |6|. We focus on this regime be- 
cause the wavefunction in this limit acquires an analytic 
structure which allows exact evaluation of the energy of a 
vortex lattice. As a result, it is possible to scan through 
a wide range of lattice structures which would be im- 
practical for numerical calculations because of the time 
and the accuracy required. Although not directly appli- 
cable to current experiments on vortex lattices (which 
are performed at lower angular momenta), the physics of 
the mean field quantum Hall regime is still quite relevant 
as the vortex lattices in these two regimes are connected 
continuously to each other. 



One special feature of the majority of two-component 
gases so far studied, (notably mixtures of hyperfine states 
of ^''Rb 1^ in magnetic traps, or ^'^Na |^ in optical traps) 
is that the interactions between like species (denoted gi 
and (72) and unlike species (denoted 312) are very simi- 
lar, within a few percent of each other. Thus, if there 
are an equal number of bosons in each component, and if 
the trapping potentials of the two components are made 
identical, then the two components will be the same size 
and contain the same density of vortices. In this case, 
one expects that each component will contain identical 
vortex lattices, with one lattice displaced relative to the 
other. While we are mainly interested in the experimen- 
tally relevant cases, where 51 ~ 32 ~ 512, considerable 
insight is gained by studying vortex lattices as a function 
of the interactions. Considering the case 51 ~ 52 7^ 512, 
we find a wide range of vortex lattice structures as the 
parameter a — 3x2! ^g\gi is varied. The vortex lattice 
has a fixed structure over certain intervals of a, while 
it varies continuously in others. The structure near the 
isotropic point g\ = g2 = gn consists of identical rectan- 
gular lattices in both components, with one displaced to 
the center of the unit cell of the other. The aspect ratio 
of the unit cell changes with a, and is \/3 when a ~ 1. 

The mean-field quantum Hall regime: The con- 
densate wavcfunctions and '52 of a two-component ro- 
tating Bose gas are determined by minimizing the grand 
potential K = E — VlL^ — fiiNi — H2N2, where E is the 
energy of the system, fl is the rotational frequency, Lz is 
the angular momentum along z, and fj,i (i = 1, 2) are the 
chemical potentials fixing the number of bosons A^i and 
N2 in each component. For simplicity, we assume identi- 
cal trapping potentials for each component. We consider 
a cigar shaped trap with the symmetry axis z coinciding 
with the axis of rotation. As discuss in Q , the slow varia- 
tion of the trapping potential along z allows one to apply 
a Thomas-Fermi approximation for the z dependence of 
'i'i and write K a.s J dzdr /C(r, z). 
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with r = {x,y), Hi{z) = - ^MujIz^, = Anh'^ai/M, 
{i = 1,2), and 512 = 47r?i^ai2/M, where a^, and aij 
are the s-wave scattering lengths between hke and un- 
hke bosons respectively. As z is treated as a param- 
eter, it is convenient to write '^i — y/ni{z)^i{r; z), 
with J \^i{r; z)\^ (fr = 1. The number constraint 
/ drdz l^fp = Ni becomes J ni{z) dz = Ni. 

As il approaches Wj^, the wavefunctions <i>i are made 
up of the orbitals Wm(r) in the lowest Landau level in the 
xy-plane, $i(r, z) = J2m=o ^rn{z)um{r) , where Ujnjr) = 
(27rm!)-i/2[(x + i?/)/d]"e-'''/2d^, and d = ^Ti/Muj^. 
The potential K, then becomes 
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where (r^)i = /r^|$ip o?-^r, and 
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As shown in ref. wavefunctions in the lowest Landau 
level (not normalized) can be written as 



(r) = A - aa)e 



ly. 
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where A is an arbitrary constant and {a^} are the zeros 
of (f). If the zeros form a infinite lattice with unit cell size 
Vc, it is shown in Q that |</)p is a product of a Gaussian 
and a function periodic under lattice translation, i.e. 



5(r), g(r)=5(r + R), 



(7) 



where R = n-iBi +71262, ni are integers, and Bi, B2 are 
the basis vectors of the lattice. The width a reflects the 
number of vortices of the system. It is given by 



(8) 



The periodicity of g{r) imphes g{v) — X^k^k^'^ '^j 
where {K} are the reciprocal lattice vectors. 

In the following, we shall consider a two-component 
Bose gas with equal particle numbers and trapping po- 
tentials, and with interactions gi ^ g2 ^ <7i2- If 5i = <?2, 
the two components are identical and one expects that 
each will contain identical vortex lattices, translated with 
respect to one another. Sufficiently small differences in 
<?i ~ 92 should not change this structure (though changes 
may occur in the density profiles rti(z), the parameters 
of the lattice, and the relative displacement rg). This 
structure persists because, even when gi 7^ 52, the two 
components contain equal vorticity, hence an equal den- 
sity of vortices. The potential energy is minimized by 
interlacing the two lattices; if the vortex lattice in one 
component were to deform, the other has to follow to keep 
the interaction energy at a minimum. We shall therefore 



consider normalized condensates $1 and 4>2 in Eq- @) 
with densities of the form 

l*iP = ('ra')"'EK5Ke^''-""e-'^'/"' (9) 
|$2p - (7ra2)-i EK5Ke^^-('--°)e~'-V<^^ (10) 

5K = WfEK'5K'e-'^'K"/4). (11) 



The wavefunction is described by variational parameters 
ni{z), tj^, the basis vectors B^ (which determine the unit 
cell size Wc), and the relative displacement Vo- 

By integrating Eqs. (||) and (^, one sees that up 
to terms of relative order Wc/c^ the cloud's radius is 
(r^)i — {r^)2 — cr^- Defining the quantities I and /12 as 
jm^d^r = I/ina^) and /|<I>i|2|cI>2|2 d^r = /i2/(7ra2), 
we have 
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h2 = 2^.gK5K'e °e ' ^ ' ' , 



(12) 
(13) 



K 



and the potential /C takes the form 

JC= - [fi{z) - hui^ - h{uj^ - n){(7^/d'^)] (ni + 712) 

+ (c?V2ct2) (^njgil + nlg2l + 2nin25i2/i2) • (14) 

To evaluate /C, we need an explicit expression for the 
coefficients ■ These coefficients are most easily derived 
by introducing the complex representation for the basis 
vectors, hi = (x -f iy) ■ B^. The area of the unit cell is 
then = i{b\h2 — h2hi) /2. If we orient the lattice so that 
hi is real, i.e. Bi = 6ix, B2 = 5i(ux + fy), we then have 



h2 = hi{u + iv), v^ = h\v. 



(15) 



In the Appendix, we show that a Junction $ in the low- 
est Landau level describing a regular vortex lattice con- 
tained in a cylindrically symmetric cloud will have the 
form (j){T) — f{w)e~^ 1'^'^ , with w = x -\- iy, and 



f{w)^9iC,T)e^^'/^^^ 



(16) 



where C. = w/bi = {x -\- iy)/bi, t = u-\- iv = b^/bi, and 
is the Jacobi theta- Junction defined as 

-. 00 

6{C,t) = - y (_l)»e^-T(»+l/2)%2.^C(n+l/2)^ (^7) 



The density \(f>^ is therefore of the form Eq. (j^, with a 
given in Eq. M), and 



?(r) = \diC,T)expi-7ry^/v,)\' 



(18) 



The Jacobi theta-function has the quasi-periodic proper- 
ties 



d{c + i,T) = e{c,T), 

^(C + r,r) = -e— (-+2«)0(C,r), 
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which imphes the periodic property g{r) = g(r + R). The 
Fourier coefficients of ^(r) are 



<?K = (-1) 



mi +m2+mim2 — t'^ |Kp /Stt 



(21) 



where K = miKi + 7112^2, and are the basis vector 
of the reciprocal lattice, Ki = {2tt/vc)'B2 x z, K2 — 
{2tt/vc)z X Bi, and 



(2^)2 



((wmi)^ + (to2 — wmi)^) 



(22) 



Since we work in the limit of large vortex number, the 
size of the cloud is much larger than the unit cell, i.e. 
TTa'^ /vc >> 1. We can therefore ignore all K + K' 7^ 
terms in Eqs. 
then have 



(|T|) and (|l|), since a^K^ > na^/vc. We 
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where gj^ is given by Eq. ( pi]) and K-sum is over the 
integers mi,m2. Since the expressions of / and I12 in 
Eq. ( p3|) are independent of cr^, the minimization of K. in 
Eq. ( ]l4| ) with respect to cr^ and become very simple. 
The optimum ct^, n = (ni, ri2), and /C are given by 



= cP iniz) - huj^_) I {3h{uj^ - O)) , 
n(z) = {2/3){a^/d^)(p{z) - hu^) ■ 1, 
K = -{l|Z){^l{z) - huj^) 1 • n(z), 
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It is clear from Eqs. ( p^ through (^^ that the solution 
for the case where gi — 92 « \gi +92] is very close to 
that of gi =52- The lattice shape (parameterized by 
To, u, and v) enters the grand potential only through the 
factor 1 • G^-*" ■ 1. When gi — g2 this term is inversely 
proportional to J = / + a/12, and the most favorable 
lattice is the one that minimizes J. 

Summary of Results: It is interesting to compare 
the two-component case with the single-component case. 
In the latter system, energy minimization reduces to min- 
imizing /. The only local minimum is the triangular lat- 
tice, where / = 1.1596; the square lattice is a saddle 
point with / — 1.1803. The minute difference between 
these values of / makes a simple numerical minimization 
of (|l|) challenging and illustrates the utility of the ana- 
lytic scheme used here. 

For a two-component Bose gas, the most favorable lat- 
tice minimizes / -I- a/12 ■ In the minimization it is conve- 
nient to measure lengths in unit of the the basis vector 
Bi = 61X, and write complex representation of B2 and 
To StS T = u + iv = |r|e*'' and To = a + br respectively. 
The phase diagram of the vortex lattice as a function of 
the ratio a — 312/3 is shown in Figs. 1 and 2. The major 
features are 



(a) a < 0: In this region the vortices of the two com- 
ponents coincide with each other (a = 6 = 0) to form a 



triangular lattice (r 



(b) < a < 0.172: the vortex lattice in each component 
remains triangular. However, Tq undergoes a first order 
change so that one lattice is displaced to the center of 
the triangle of the other (a = 6 = 1/3). The lattice type 
(characterized by t = e™'^) remains constant within this 
interval. 

(c) 0.172 < a = 0.373: At a = 0.172, jumps from 
the center of the triangle (i.e. half of the unit cell) to 
the center of the (rhombic) unit cell {a = b = 1/2). The 
angle 77 jumps from 60° to 67.95° at a = 0.172, and 
increases continuously to 90° as a increases to 0.372. As 
a result, the lattice type is no longer fixed and the unit 
cell is a rhombus. The modulus of r, however, remains 
fixed across this region. 

(d) 0.373 < a < 0.926: The two lattices are "mode- 
locked" into a centered square structure throughout the 
entire interval {t ^ i, a = b ~ 1/2). 

(e) 0.926 < a: The lattice type again varies continu- 
ously with interaction a. Each component's vortex lat- 
tice has a rectangular unit cell (77 = 7r/2) whose aspect 
ratio |r| increases with a. Both ®^Rb and ^^Na have 
interaction parameters within this range. At a = 1, 
(51 = 92 = gi2 = g), the aspect ratio is -\/3- If one 
ignores the difference between the components, the com- 
bined lattice is triangular, as is expected. 

It is interesting to note that in the absence of rotation, 
the two components change from miscible to immiscible 
when a increase beyond 1. No such change, however, 
happens at a = 1 in the high angular momentum limit. 
This qualitative difference in behavior occurs because the 
presence of a vortex lattice naturally modulates the den- 
sity of each component, with the high density regions of 
one fluid coincident with the low density regions of the 
other. Thus the system is effectively phase separated 
whenever staggered vortex lattices are present, even for 
a < 1. In particular, the vortex lattice near a = 1 (above 
or below) is made up of alternating rows of vortices of 
each component, (see Fig. 1), and the system therefore 
contains stripes in which one component has a high den- 
sity and the other component has a very low density. As 
a increases, the stripes become more pronounced. 

Final Remarks: The diversity of the vortex lattice 
structures in the two-component Bose gas has once again 
demonstrated the rich properties of these systems. Our 
calculation, based on exact evaluation of the vortex en- 
ergy, assumes a perfect lattice. Considering the long re- 
laxation times in clouds of dilute atoms, one might see 
more complicated structures, where patches of vortex do- 
mains are separated by defects or grain boundaries. Nev- 
ertheless, the underlying equilibrium structure should be 
reflected within each vortex domain. 

So far, we have only discussed two-components sys- 
tems with simple interpenetrating Bravais lattices. Our 
method is more general in that it allows the exact eval- 
uation of the energy of an arbitrary regular lattice (with 
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however complicated unit cell decoration). Such struc- 
tures may be favored when the particles in each com- 
ponent have different numbers, trapping potentials, or 
masses (as in the case of ^^Na-^^Rb mixtures). 

Appendix: The general form of a vortex lattice in 
the lowest Landau level is '9{x,y) = f{w)e~'^ 1'^'^ , where 
w = X -\- iy and / is an entire function whose zeros form 
a regular lattice {6 = nh\ ~\- 71-2^2}, where Ui are integers, 
and 61 and 62 are the complex basis vectors, (&2 — h\(u~\- 
iv)). Since the Jacobi theta function 0{{x + iy)/bi, u + iv) 
is an entire function with exactly these zeros, we have 
/(w) = e{C,T)h{C), where ( = (x + iy) /hi = x + iy, 
T = u+iv = 62/&1, and h{(^) is an entire function without 
zeros. To ensure the normalizability of ^P, this function 
can only be of the form h{C,) = exp(ciC + C2C^)- It is 
straightforward to show that 

l^'(C, t)\^ = ^„^(-l)"'e2'^""*e-'^""'/2i„ (28) 

j-^ jj_ ^ ^ ^-|^ gZ7r{m+7n')^ ^(iTrum — 27ry — '7ru7n'/2)m' (29) 

By applying the Poisson summation formula to to Eq. 
(P9|), we have 

= y?^(-i)(™+i)'=e(-"('=+"'"+''^)'/2.). (30) 
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FIG. 1. Phases of the two-component lattice: black and 
grey dots represent vortices of each of the two fluids. The 
panels a) through e) show the vortex structure in each of 
the phases described in the text. The final panel depicts 
the geometry of the lattices; the black and grey dots re- 
spectively occupy positions in the complex plane {m + nr}, 
and {(o -I- m) -I- (& -I- ?i)t}, where m, n are integers. All mini- 
mal-energy configurations have a — b. 



We thus have 



g27rj/^/«c 



(31) 



where r = xx + yy, K = (2Trmx — 2TT{n+um)/vy)/bi, and 
(7k is given by (pi|). The density of the system is then 
|5-(r)|2 = |6l(C,r)|2|e'=i';+'=^«'|2e-'''/'^'. For a vortex lat- 
tice with inversion symmetry about the origin r = 0, we 
have ci = 0. In addition, if the cloud's envelope is cylin- 
drically symmetric, we have C2 = 7r/(2uc), which gives 
Eqs. (0), (H), ( p^ and (jl^). Similar approaches have 
been used by Tkachenko ||9|| and Abrikosov in their 
respective studies of ''He and type two superconductors. 
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FIG. 2. The parameters of the vortex lattice as a function 
of a = gi2 1 yjgigi, a measure of the importance of interactions 
between unlike atoms. The phases, labeled a through e are 
illustrated in Fig. ^ along with the parameters t = |T|e"' 
and a. Solid and dashed vertical lines respectively denote first 
and second order phase transitions. The open circle on the 
horizontal axis indicates a = 1. 
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